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The Open-System Dicke-Model Quantum Phase Transition with a Sub-Ohmic Bath 
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We show that the critical exponent of a quantnm phase transition in a damped-driven open system 
is determined by the spectral density function of the reservoir. We consider the open-system variant 
of the Dicke model, where the driven boson mode and also the large N-spin couple to independent 
reservoirs at zero temperature. The critical exponent, which is 1 if there is no spin-bath coupling, 
decreases below 1 when the spin couples to a sub-Ohmic reservoir. 
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Quantum critical behavior appears in driven dissipa¬ 
tive systems when the steady-state [1-3] , rather than the 
ground state of a Hamiltonian [4-6], undergoes a non- 
analytic, symmetry-breaking change at a critical param¬ 
eter value. The interplay of an external coherent excita¬ 
tion and the dissipation can lead to a steady-state which 
is far from the ground or thermal state. Driven dissi¬ 
pative systems cannot be, in general, mapped onto an 
effective Hamiltonian system. It is thus unclear how the 
critical behavior of the open system is related to uni¬ 
versality classes of known quantum and thermal phase 
transitions [7]. 

Recent observation of the Dicke-model superradiant 
phase transition motivates us to raise this question. Ul¬ 
tracold atoms coupled to the radiation field of an opti¬ 
cal resonator allowed for the quantum simulation of the 
Dicke model [8-10] and the experimental demonstration 
of the phase transition [10-16]. The boson component 
of the model is represented by a single mode of a high 
finesse optical cavity. The spin-A^ component is effec¬ 
tively realized by constraining the motion of ultracold 
atoms into the space of two momentum eigenstates. The 
interaction is implemented by a far-detuned laser field il¬ 
luminating the atoms from a direction perpendicular to 
the cavity axis. Photons scatter into the cavity, which 
has a recoil on the atoms. Above a threshold of the laser 
intensity, which translates to the coupling strength be¬ 
tween the spin and the boson mode, a mean field of the 
cavity mode and the large spin is formed spontaneously. 

Since the cavity mode is coupled through the mirrors 
to the outside electromagnetic vacuum field, this system 
is intrinsically open. It is also a substantial feature that 
the coupling between cavity photons and atoms is gener¬ 
ated by an external laser. In the frame rapidly rotating 
at the laser frequency, the time-dependent driving can be 
eliminated and the remaining low-frequency dynamics is 
described by the time-independent Dicke-type Hamilto¬ 
nian [9, 10, 17]. Note, however, that the outcoupled field 
is continuously supplied by the external laser and a pho¬ 
ton current is driven through the system. 

According to the effective Hamiltonian, the incoherent 
photon population in the ground state diverges at the 


critical point as a power law with exponent 1/2 [18]. At 
variance, when the cavity loss is taken into account [19, 
20], the exponent changes to 1 which is much closer to 
the experimentally measured value of about 0.9 [21]. The 
dissipation and the accompanying quantum fluctuations, 
even at zero temperature, has thus been demonstrated 
to substantially modify the correlation functions and the 
critical exponents. 

We will show in this Letter that the exponent can be 
continuously tuned below 1 and its actual value is de¬ 
termined by the low-frequency spectral properties of the 
bath. The value 1 follows from the Markovian dynamics 
imposed by photon loss: the ‘real’ frequency of the leaky 
photons is in the optical range where the electromagnetic 
field has practically a constant spectral density function 
in the relevant narrow, say, few kHz bandwidth. The 
effect of such a reservoir on the coupled system can be 
interpreted effectively as that of a thermal bath at a non¬ 
zero temperature on a single boson mode [22]. However, 
the bath can have a non-trivial spectral density function. 
For example, in the above-quoted experimental realiza¬ 
tion there is another dissipation channel acting on the 
spin component of the system, which was shown to have 
a significant influence on the measured correlation func¬ 
tions [21]. Based on a microscopic calculation beyond the 
Bogoliubov approach, the observed damping has been at¬ 
tributed to a Beliaev-type damping process in the super¬ 
fluid with a highly non-trivial spectrum of the phonon 
reservoir at low frequencies [23, 24]. 

We consider a bosonic two-mode model in which one of 
the modes is coupled to a simple Markovian bath whereas 
the other one is subjected to a colored reservoir. This 
general model describes with high accuracy the normal 
phase of the Dicke-type phase transition [9] . The generic 
feature is that the soft mode of the phase transition is 
composed of two bosonic modes which are coupled by 
the Hamiltonian 

Hs/h = Wao)a + ujbh^b + |(a -I- a'^)(& + b^) ■ (1) 

The coupling y is restricted to the range below the criti¬ 
cal point yc, which we will determine later. The Keldysh 
path integral approach is invoked to calculate the dis- 
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sipative effects [22, 25]. The action of the Markovian 
oscillator in frequency space reads 
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where Oci(w) and aq(w) are the classical and quantum 
fields corresponding to the mode a. This mode a is like 
the ‘photon’ mode in the Dicke model, has a high fre¬ 
quency and iOa is referenced to some driving frequency. 
This mode emits then “high frequency” photons into the 
vacuum, thus its decay is unaffected by the interaction 
between modes a and b. Moreover, the flat reservoir 
spectrum at high frequencies ensures the validity of a 
Markovian approximation, which is reflected by using a 
single constant parameter k. It is important to note that 
the time dependent external driving, formally eliminated 
from the Hamiltonian (1) in a procedure presented in 
detail in Ref. [17], has as a consequence that a Marko¬ 
vian decay with frequency-independent self-energy func¬ 
tion can be used in the above action. 

The genuine low-frequency mode is b, which interacts 
with a sub-Ohmic bath at zero temperature. It is essen¬ 
tial that this reservoir is sampled at the eigenfrequencies 
of the coupled system. The dissipative action that de¬ 
scribes the dynamics reads 



( 3 ) 


where 6 ci(w) and 5q(w) are the classical and quantum 
fields corresponding to mode b. The reservoir is char¬ 
acterized by the coupling-density profile p{uj) which de¬ 
fines, at zero temperature, the Keldysh component of the 
action, D{u!) = 2iTrp(uj), and also the retarded and ad¬ 
vanced level-shift functions. 


K^^^(oj)=V [ ——^ do; ^iTrp{uj). (4) 
Jq lu — lu 

Note the symmetry K^{uj) = [K^(u!)]*. For a sub- 
Ohmic reservoir, p{uj) starts with an exponent 0 < s < 1, 
namely 


p{u}) = 0 (w)- 
TT 
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( 5 ) 


where 7 is the dissipation strength, ojm is a cutoff fre¬ 
quency, and ©(w) is the Heavyside function. This reser¬ 
voir yields a non-conventional level-shift function, which 
is purely real for w < 0 and complex for w > 0. The 
analytical continuation of K^(pj) to the upper half, and 
K^{uj) to the lower half of the complex plane can be 
fitted together in the function 


K{z) 
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such that it has a branch cut and renders = 

lim^_j.o ± irf} on the positive real axis. 

The Keldysh action that corresponds to the interaction 
term in Eq. (1) is 

/ [(<■.+“;)('>=■+‘a++Of*'.+'>;>]. 

( 7 ) 

which contains counter-rotating terms, like Udbq 
etc. Following the method of Ref. [22] , the variable space 
is doubled by introducing fields with negative frequencies 
so that the total Keldysh action S = Sa + Sb + Sab can 
be expressed in a closed quadratic form 



with the 8 -component field 


v'''(w) = [o*;(a;),aci(-a;),5*;(a;),6c/(-w), 

a*q{u}),aq{-u}),b*{uj),bq{-uj)]. (9) 

The 4x4 blocks are matrix Green’s functions. The char¬ 
acteristic frequencies of the system correspond to the 
poles of the retarded Green’s function continued ana¬ 
lytically to the lower half of the complex plane. The 
extension is not unique, we use the level-shift function in 
Eq. ( 6 ) on the second Riemann sheet 

-ITTS 
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and require the symmetry property K^{z) = [K-^{z)Y. 
The poles come then either in pairs (zi, —Zq ), i.e., op¬ 
posite real part and the same imaginary part, or they 
are purely imaginary. The poles obey the characteristic 
equation det[GK^ 4 ]“^(z) = 0 , i.e.. 


[(z + z«)2-a;^][(z-zF(z))2-(a., + A(z))2] 

- -I-A(z)) = 0, (11) 

where r(z) = {K^{z) — KTff(—z*))/( 2 z) and A(z) = 
(iFif(z) + i^fl(-z*))/ 2 . 

The numerical solution of this equation for the pole 
corresponding to the soft mode frequency uig is shown 
in Fig. 1 as a function of the control parameter y for 
various values of the reservoir coupling strength 7 . For 
ujb < UJa, the soft mode grows, as y is increased from 
0, out of the bare mode b. A characteristic feature of 
second order dissipative phase transitions is that the real 
part of the soft mode frequency (top panel) decreases 
to zero first, and at this exceptional point a linewidth 
bifurcation takes place [26]. The larger the 7 , the real 
part vanishes for smaller y. The critical point is the one 
where the upper branch of the imaginary part (bottom 
panel) reaches zero. 


Vc = 
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FIG. 1: (Color online) Real part (upper box) and imaginary 
part (lower box) of the soft mode frequency uis as a function of 
the coupling y to the sub-Ohmic bath for various dissipation 
strengths, 7 = 0 (solid red), 0.1 (dashed green) and 0.5 iJb 
(dotted blue). The bath low-frequency exponent is s = 4/5. 
Other parameters are LUa = k = ‘iuib- 

regardless the coupling strength 7 to the bath. This is be¬ 
cause the level-shift function vanishes at zero frequency, 
r(z —>• 0 ) = 0 and A (0 —>• 0 ) = 0 , hence it drops out 
from the z = 0 solution of Eq. (11). The vanishing of the 
soft mode frequency implies that there is a divergence in 
the correlation functions at y = yc- The exponent, how¬ 
ever, cannot be extracted merely from the position of the 
poles, because the dynamics in a coloured reservoir is not 
Markovian. Therefore, we resort to a direct, numerically 
exact calculation of the correlation functions. 

The Fourier-transform of the correlation functions, 

C,(t) = ({a(t),at(0)}), a(t) = ({5(t),6t(0)}), (13) 

are given by the components of the Keldysh Green’s func¬ 
tion, Ca{u}) = i[Gf^ 4 ]ii and Cf,(w) = respec¬ 

tively. It can be obtained as 

Gfx4 = -Gfx4(^)Dfx4(^)Gfx4(w) . (14) 

Figure 2 presents the correlation functions Ca{uj) (a) 
and Cb{uj) (b) for different values of the coupling y ap¬ 
proaching j/c- The former one is measured directly in the 
ultracold atom realization of the Dicke model, since it is 
the power spectrum of the field leaving the cavity. For 
a non-interacting system {y = 0 ), the photonic correla¬ 
tion function (Fig. 2a, solid red) is simply a Lorentzian 
peak around uj = uja oi width 2k. It has a non-vanishing 
part for negative frequencies, causing effectively a heat¬ 
ing, which is due to the underlying driving at optical 
frequency. The correlation function of mode b for y = 0 
(Fig. 2b, solid red) reflects the properties of the colored 
bath. This peak is strictly zero for w < 0 and has a tail 
only for large uj, which is consistent with the choice of 
the coupling density profile, Eq.(5). Though the spec¬ 
trum is not Lorentzian, the resonance peak corresponds 
well to the real part of the soft-mode frequency, which is 
significantly shifted from the bare frequency Wh = 1 to 
ojh = 0.5 in accordance with Fig. 1 (dotted blue). 
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FIG. 2: (Color online) Power spectrum of the ‘photonic’ mode 
a (a) and of mode b (b) for various coupling strengths ap¬ 
proaching the critical point: y = 0 (solid red), 0.5 (dashed 
green), 0.9j/c (dotted blue) and 0.95yc (dashed-dotted brown). 
Parameters: uJa = k = 2oJb, 7 = 0.5 ca;,, s = 4/5. In the pho¬ 
tonic spectrum (a), the peak for y = 0.95j/c is divided by a 
factor of 10 to fit into the plotted range. 

On increasing the coupling y, the resonance peak of 
mode b moves towards w = 0, following the behaviour of 
the soft-mode in Fig. 1. The dotted blue line presents 
the correlation function slightly above the point where 
the real part of the soft mode vanishes {y « 0.93?/c)- 
Such close to the critical point, the correlation function 
cannot be simply understood in terms of the soft mode. 
This point y « 0.932/c demarcates two different kinds of 
correlation function also for the mode a, shown in panel 
(a). Below, there is a doublet symmetrically to w = 0 
on top of the broad Lorentzian corresponding to the soft 
mode frequency shown in Fig. 1 (the opposite sign eigen- 
frequency is not shown there). This doublet uncovers 
the hybridisation of mode a with mode b. For a larger 
coupling, e.g. y = 0.95j/c (dashed-dotted brown), a sin¬ 
gle sharp peak develops in the origin, that diverges for 

y Vc- 

The mean steady-state population (a^(O)a(O)) = 

{C{t — 0) — l)/2, i.e., cavity photon number, is given 
by the integral 

Ca{t = 0) = J ^Caito) . (15) 

The photon number diverges at the critical point with a 
scaling exponent, which is determined numerically. The 
inset of Fig. 3 shows the integrated correlation function 
Eq. (15) for variable coupling 1 — y/yc on a double log¬ 
arithmic scale. In the close vicinity of the critical point, 
\y ~ J/cl < 10 “^, the power-law can be recognized and 
the exponent can be extracted. The main result of this 
Letter is summarized in the main panel which presents 
the continuous and monotonous variation of the critical 
exponent as a function of s. For a sub-Ohmic bath the 
exponent decreases below 1 . Note that the ohmic bath 
with s = 1 renders the critical exponent 1 which was 
known also for the Markovian dissipation without the 
colored bath (7 = 0) [19, 22]. 

The crucial role of the sub-Ohmic bath in the criti- 
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FIG. 3: (Color online) Critical exponent as a function of the 
exponent s of the sub-Ohmic bath, which is extracted from 
the scaling of the population in mode a (aka ‘photon number’, 
solid red line) and in mode b (green circles). The power-law 
scaling of the photon number is shown in the inset for s = 0.8 
(red crosses -f), 0.7 (green crosses x), 0.6 (blue stars), 0.5 
(brown squares). Other parameters are the same as in Fig. 2. 
The solid black line shows the 7 = 0 case for reference. 

cality of dissipative systems appears also, for example, 
in the spin-boson model [27], where the exponent of 
the localization-delocalization transition depends on the 
low-frequency exponent of the spectral function [28-30]. 
The transient (out-of-equilibrium) dynamics exhibits a 
crossover from the delocalised to the localised fix points 
of the equilibrium phase transition [31, 32]. The Ising 
model criticality is also non-trivially influenced by cou¬ 
pling the spins to a bosonic bath [33]. In these cases, the 
ground state of the coupled system is in question. By 
contrast, here we considered an externally driven system 
in which the dissipation, more precisely, the interaction 
with a bath, leads to a steady state. There is a flow of en¬ 
ergy through the coupled a and b bosons subsystem from 
the coherent driving into the reservoir, which cannot be 
accounted for by an effective Hamiltonian. The critical 
point in the steady-state corresponds then to a genuinely 
non-equilibrium quantum phase transition [34, 35]. This 
Letter was the first to study non-Markovian bath effects 
in a non-equilibrium quantum criticality. We calculated 
in a numerically exact way the critical exponent in the 
open-system Dicke phase transition, and showed that it 
is a monotonous function of the low-frequency exponent 
of the bath spectrum. 
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